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O . Abstract 



The alleged mistakes in recent papers that reanalyze the backgrounds to the 'LSND 
anomaly' do not exist. We maintain our conclusion that the significance of the 'LSND 
anomaly' is not 3.8 a but not larger than 2.3 a. 
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In two recent papers (HIS, we reanalyzed the backgrounds to LSND's signal of — > u c 
oscillations (the 'LSND anomaly'), questioned LSND's analysis procedures, and claimed 
that the significance of the 'LSND anomaly' is not 3.8 a but not larger than 2.3 a. 

The concerns discussed in our papers are: 

1. inadequate knowledge of n ± production by 800 MeV/c protons on various target nuclei 
at the time when LSND calculated their neutrino fluxes; 

2. inadequate or missing consideration of pion production by higher-generation protons, 
pions, and — in particular — neutrons; 

3. the inclusion of the first bin in the fit of the R y distribution that does not carry any infor- 
mation on the distinction between correlated and accidental 7's; 

4. questions on the 'efficiencies' of correlated and accidental 7's; 

5. questions on the effective rate of accidental 7's and their influence on the i? 7 distributions 
of correlated and accidental 7's; 

6. missing systematic errors of the 'base distributions' of correlated and accidental 7's; and 

7. missing positrons from 12 N gs beta decays that are misidentified as correlated 7's. 

Of these seven items, only items 1 and 7 are addressed in the paper 'Corrections to the 
HARP-CDP Analysis of the LSND Neutrino Oscillation Backgrounds' by G.T. Garvey et al. [ffl 
that claims several mistakes in our analysis. 

Before we argue in the following that these alleged mistakes do not exist, a clarification 
is in order. LSND insist on the notion that their neutrino flux calculations are correct within 
their quoted errors. We disagree with this notion. First, the knowledge of pion production 
cross-sections has considerably improved in the last decade (we use state-of-the-art Geant4 and 
FLUKA cross-sections, further improved by pertinent experimental results that are reported in 
Ref. [1]). Second, LSND apply the same percentage errors to neutrino flux integrals and to parts 
of the neutrino spectra — in particular to the small high-energy portions of the spectra, which is 
not realistic. Third, the LSND parametrization of pion cross-sections does not describe the 
results of their calibration experiment E866 0. All this leads to much too optimistic error as- 
signments on neutrino fluxes by LSND, and to misconceptions on what quantities have actually 
been 'measured'. 

Our simulation of LSND's setup serves the following purpose: to assess quantitatively the 
effect of the improvement of pion production cross-sections over those that were employed 
by LSND more than a decade ago. This necessitates (i) a simulation that can switch between 
the application of the LSND cross-sections ('LSND emulation') and our 'best estimate' cross- 
sections; and (ii) a demonstration that this simulation reproduces, with LSND cross-sections, 
approximately the published LSND result on the conventional v c background. 

Our 'LSND emulation' reproduces with 0.585 x 10~ 12 (pot.cm 2 )" 1 remarkably well LSND's 
published result on the conventional v c background of 0.65 x 10~ 12 (pot.cm 2 )" 1 . Then we 
add the observed increase when switching from the 'LSND emulation' cross-sections to our 
'best estimate' cross-sections, to the conventional u c background that LSND published. In this 
procedure small differences between the real LSND setup and its simulation in our program 
cancel, and we obtain the conventional u c background that LSND should have obtained with 
our 'best estimate' cross-sections. The stringent result is an increase by a factor of 1.6. 

The u c /u^ ratio is not a useful quantity for the normalization of neutrino fluxes as u c orig- 
inate from pT from n~ decays while originate from /i + from 7r + decays. In the LSND 
energy domain n + production is essentially unrelated with n~ production. Therefore, it makes 
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no sense to introduce an unnecessary error by transporting the uncertainty in ir + production into 
the calculation of the conventional u c background which originates from 7r~ only. 

Our 'best estimate' of and rates from pion decay in flight is larger than the LSND 
estimate by factors of 2.5 and 3.3, respectively. The and rates refer to E v > 123.7 MeV for 
Up and to E v > 113.1 MeV for 9^ and hence to very small high-energy portions of the respective 
neutrino spectra. These portions are not relevant for the calculation of the conventional z7 e 
background. Rather, they are relevant for (i) calculating what we termed 'Background IF that 
concerns primarily events where a muon is misidentified as electron, and (ii) to appreciating 
LSND's constraints on neutrino fluxes (which in no way contradict our claim of a conventional 
u e rate larger by a factor of 1.6). The factor of 3.3 for is essentially unrelated with our 
claim of a v e background larger by a factor of 1.6 than calculated by LSND. The proportionality 
between these two quantities purported by Garvey et al JU lacks justification. 

We agree that '3 MeV is not a hard cutoff for muons that are not observed. Rather, it is 
a reasonable estimate. It is for this reason that '3 MeV is nowhere used in our calculation 
of 'Background II' . We consider it safer to calculate Background II from the small portion of 
the neutrino spectrum between the neutrino energy threshold and 4 MeV above. The value of 
4 MeV is approximate and beset with instrumental uncertainties; we bypass these uncertainties 
by considering not absolute numbers but relative numbers: we determine a factor by dividing 
'best estimate' predictions by 'LSND emulation' predictions, and apply this factor as a correc- 
tion to LSND's estimate of Background II. 

The final LSND physics paper @ comprises a discussion of backgrounds that mimic 2.2 
MeV 7's. There, even backgrounds as low as 0. 1 events are discussed. There is no mention of a 
background from 12 N gs beta decays that is now claimed to be ~0.2 events. As for the claim that 
an electron produces many more PMT hits than a photon, we note that the energy of a 2.2 MeV 
photon is in mineral oil within a few cm fully converted into the kinetic energy of electrons 
by Compton scattering and photoelectric interaction. The R distribution of betas from 12 N gs 
decays can neither agree with the i? 7 distribution of correlated 7's nor with the FL, distribution 
of uncorrected 7's, for the pulseheight distribution of the beta's is intrinsically different. It is 
straightforward to calculate the background from 12 N gs beta decays that mimic correlated 7's, 
and we stand by our estimate of a background of 2.3 events that was neglected by LSND. 

In summary, we maintain our conclusion that the significance of the 'LSND anomaly' is not 
3.8 a but not larger than 2.3 a. 
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